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CN ■ Abstract 

In non-equilibrium states, currents are produced by irreversible processes that take a system 

VO ' toward the equilibrium state, where the current vanishes. We demonstrate, in a general setting, 

(N' 

that a superstatistics arises when the system relaxes to a (stationary) quasi-equilibrium state 

(—1 '. 

Q . instead, where only the mean current vanishes because of fluctuations. In particular, we show 

(p : 

that a current with Gaussian white noise takes the system to a unique class of quasi-equilibrium 
states, where the superstatistics coincides with Tsallis escort ^-distributions. Considering the fully 

■^ ■ developed turbulence as an example of such quasi-equilibrium states, we analytically deduce the 

power-law spectrum of the velocity structure functions, yielding a correction to the log-normal 
__i ^ model which removes its shortcomings with regard to the decreasing higher order moments and 

Q ' the Novikov inequality, and obtain exponents that agree well with the experimental data. 
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INTRODUCTION 

Recent experimental and numerical studies have shown that many complex systems, and 
of particular relevance to this work, the fully developed turbulent system (iKZl], exhibit 
non-Gibbsian distributions that are well fitted by Tsallis g-distributions |8|. When these 
distributions are involved, physical mean values correspond to the g-average values, which are 
averages over the postulated escort g-distributions. The physical interpretation of the escort 
distributions is, however, not clear. Addressing from a physical foundation the question of 



why such generalized distributions arise is crucial 



m 



complex systems. Despite several attempts |4j.l9l-lll|. a comprehensive answer to this question 



to understanding the behavior of these 



is still lacking. In this report, we approach the problem from a general perspective that 
provides a new insight and consider, as an application, the example of fully developed 
turbulence. 

We consider a macroscopic system with a relevant extensive variable X and a conjugated 
(entropic) intensive variable P = dxS, where S{X, . . .) is the system's entropy. Let Xq, Pq, 
Sq represent the same quantities for the system's environment, which acts as a reservoir. 
The closed composite system has entropy Stot = S + Sq, and the closure condition requires 
that X + Xq = Xtot (const.). If X and Xq are otherwise unconstrained, their equilibrium 



values X and Xtot ~ ^ are obtained by the vanishing of the affinity 12| (the thermodynamic 
force) F = dxStot = P ~ Pq- Thus, if F is zero the system is in equilibrium, but if F is 
nonzero an irreversible process occurs, producing a current X that takes the system toward 
the equilibrium state (dot denotes time derivative). Since the current vanishes if the affinity 
vanishes, we can expand the current in terms of the affinity with no constant term. Assuming 
sufficiently small affinities throughout, we have X = iF, where £ > is called the Onsager 
kinetic coefficient. This relationship between the current and the affinity characterizes the 
relaxation flow of the system. Thus, since 

F = (dxP) iF ^ dxP\x iF (1) 

the equilibrium state, F = 0, is an attractor provided dxP\x = ^x^lx < 0. This is, of 
course, the stability criterion which implies that the response function (like the heat capacity, 
the compressibility, and so on) is positive. 

In the equilibrium state, the probability distribution function (p.d.f.) of the intensive 
variable, therefore, is Po(-P) = ^{P — Pq)- The extensive variable is, however, distributed 



about the mean equilibrium value X with a p.d.f. PoiX) oc w{X)e^^^^^ {ks = 1), where 
wlX) is the 'density of states'. This standard result, which forms the basis of of equilibrium 
statistical mechanics, is obtained by applying the principle of equal probability of microstates 
to the closed composite system and imposing the equilibrium condition P = Pq. The latter 
condition marginalizes the P- values, implying that the equilibrium distribution, pq{X), is 
in fact a marginal p.d.f. If, instead of the equilibrium state, the current takes the system 
to a different attractor where the intensive variable is distributed about Pq, we will have 
a stationary (a quasi- or meta-equilibrium) state whose P-p.d.f., p{P), must be used to 
marginalize the P-values. In such quasi-equilibrium states, the intensive variable fluctuates 
about the mean equilibrium value so that in these states only the mean current is zero. This 
suggests that a fluctuating current (due to the underlying noise generating mechanisms) 
drives the system to such quasi-equilibrium states, where it attains zero mean value. For 
simplicity we shall consider systems with a constant density of states so that Po{X) = 
PqC"'^""^ and {X)q = X = I/Pq. (The formulation can be easily extended to the more 
general case of w{X) oc X"~^ (a > 0), where our central result in the last section, equation 
( ITO|) . is similarly shown to follow just the same). The (marginal) quasi-equilibrium X-p.d.f. 
would then be given by the Bayes rule according to 

/•oo rOD P 

p{X)= Pe-'''' p{P)dP = e-^^p(ln— )(/P (2) 

JO JO Po 

which reduces to the equilibrium X-p.d.f. Pq{X) if p{P) = po(P) = 6{P — Pq), as it should. 



Generalized distributions of the form (12]), which have been called superstatistics 



^aiiea supers 
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13|, have found many applications in diverse physical problems p, [ill, Il4l4l6l|. From our 
perspective, therefore, a superstatistics is the most natural extension of the equilibrium 
Gibbs-Boltzmann statistics to (stationary) quasi-equilibrium sates, where only the mean 
current vanishes because of fluctuations. In complex systems, we can assume that these 
fluctuations arise from the interaction of a very large number of independent microscopic 
degrees of freedom with identical stochastic properties. The central limit theorem then 
implies that the fluctuations are Gaussian distributed. Naturally assuming that they are 
temporally uncorrelated too, we have a Gaussian white noise that contains all frequencies. 
In the sequel we show that adding a Gaussian white noise to the current changes the 
attractor and yields a new attractor state with p{P) as the F-distribution, resulting in 
the Tsallis escort q-distributions as the superstatistics. In these states, the affinity and 



the current are Gaussian distributed with zero mean. In other words, we show that a 
current with Gaussian white noise fluctuations takes the system to a unique class of quasi- 
equihbrium states, where it attains the ubiquitous Gaussian distribution with zero mean 
and the superstatistics coincides with the Tsallis escort g-distributions. This explains the 
ubiquity and the physical significance of escort distributions and, hence, of the g-averages 
as quasi-equilibrium mean values. 

The noise maybe generated, in externally driven complex systems for example, by the 
variation of a control parameter. The resulting change in the attractor and the appearance 
of a new quasi-equilibrium behavior associated with the superstatistics is then an instance of 
bifurcation. Regarding turbulent flow as such, the control parameter is the Reynolds number 

nn 

and/or the boundary roughness |17l.ll8l|. the relevant extensive variable is the kinetic energy 
(per unit mass) associated with a given spatial scale, and the current is its rate of change 
per unit time due to dissipation. The current, being scale-dependent, fluctuates in a scale- 
dependent manner (intermittency) as the value of the control parameter is increased, driving 
the system to a quasi-equilibrium state where the mean dissipation vanishes at intermediate 
scales. The new behavior associated with this quasi-equilibrium state, which is identified 
with the fully developed turbulence state, is the power-law behavior of the velocity structure 
functions. We shall analytically deduce the power-law spectrum, yielding a correction to the 
log-normal model which removes its shortcomings with regard to the decreasing higher order 



moments and the Novikov inequality J19J, and obtain exponents that agree well with the 
available data. 

FLUCTUATING CURRENT: RELAXATION TO QUASI-EQUILIBRIUM STATE 

Introducing the response function Cq = — (9pX|pp > 0, we write equation ([T]) as 

F = --F (3) 

T 

where r = Co/i is the relaxation time of the system. The unique stationary solution of this 
equation is F = 0, which corresponds to the equilibrium state at t ^ r. Equation ([3]) can 
be cast as an equation for P instead. Writing -^ = 1 + -^, we have 



d P _ IF 
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Po ' Po 

F F 
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1 P 

7(p-i) (4) 



to first order in the (small) affinity. Adding the Gaussian white noise, ^(t), with the prop- 
erties 

(e(t)) = o, {m m) = 2D5it - 1') (5) 

{D > 0) to the current X = iF is equivalent to adding it to F (with D in ([5]) redefined), 
which, in turn, is equivalent to adding it to the right hand side of equation (jlj) (with D 
redefined again). We, thus, arrive at the following stochastic differential equation for P: 

This equation yield a new attractor (stationary) state, as we now show. To this end, it is 
more convenient to cast ([6]) in the form of a standard Langevin equation by introducing 
f = In -^. We have v = f{v) + ^(t) where f{v) = — ^(e^ — 1). The corresponding Fokker- 



Plank equation is [20[ dtp{v,t) + dvj{v,t) = where j{v,t) = f{v)p{v,t) — Ddvp{v,t). The 



stationary solution, p(f), at t ^ r satisfies j = 0, yielding 

p(in ^) = [r(^)]-^ (^^)^ ^^p( — ^^) (7) 

Po g - 1 g - 1 ^0 q-^Po 

where, q — 1 = Dr > 0. The P-p.d.f. p{P) = "^p(ln -p-) is, therefore, the well known 
F-distribution. Its mean is Pq and its variance, {q — 1)Pq, vanishes as D — t- (g — ?■ 1), 
i.e. as the fiuctuations disappear, recovering the p.d.f. 6{P — Pq) of the equilibrium state. 
This identifies the noise as the sole agent of driving the system toward the quasi-equilibrium 
state. Note also, by similarly writing the standard Langevin equation ([3]) (with the noise 
added, of course) as a Fokker-Plank equation, that the quasi-equilibrium affinity, and hence 
the current, are Gaussian distributed about zero, approaching the delta function as the 
fiuctuations disappear. Substituting ([7]) in ([2]), after a simple integration we obtain the 

_ g 

quasi-equihbrium superstatistics as p{X) = Pq [1 + (g — l)Po-^] '"S which coincides with 
the Tsallis escort g-distribution, p^'^\X). Of course, g = 1 + Dr > 1, where for g = 1 
the distribution reduces to the equilibrium exponential distribution. Note that the escort 
distribution itself has emerged as the physically relevant p.d.f. and is not postulated via 
the direct Tsallis q-distribution. Consequently, the quasi-equilibrium mean value {X) = 
Jl^ Xp{X)dX = (2-^)p , which exists provided g < 2 (thus, 1 < g < 2), is inherently a 
g-average. Finally, introducing the normalized variable x = X/{X), we write 

1 r /o-l\ 1"^ 
Pi^) = n 1 + I a; (8) 



and note that X may be taken to denote an extensive quantity per unit mass, when appro- 
priate. (We may note, in passing, that the physical significance of the extensivity parameter. 
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There, the authors obtain the direct 



g, coincides with interpretation given in reference 
q-distributions by considering fluctuations of the parameters of the usual exponential distri- 
bution. Their postulated Langevin equation (11), which yields the F-distribution, coincides 
with our equation (|6]), provided we make the identification x = 1/-P ^md = 1/tPq. How- 
ever, the integrand of their integral expression (4) slightly differs from ours in equation 

(EDO 

FULLY DEVELOPED TURBULENCE AS A QUASI-EQUILIBRIUM STATE 

In a turbulent flow, the current responsible for relaxation is the rate of change of the 
kinetic energy per unit mass associated with a given spatial scale. This is because of dissi- 
pation, which causes a difference between the in-fiux from the larger scales and the out-fiux 
to the smaller ones. For sufficiently high values of the control parameter (Reynolds number 
and/or boundary roughness), the system relaxes to a stationary state where the mean in- 
and out-fluxes are equal at flnite intermediate scales and dissipation occurs at a mean rate 
e only at the smallest (Kolmogorov) scale. This state, which is the state of vanishing mean 
current, corresponds to the quasi-equilibrium state of fully developed turbulence. 

Let u{r) denote the longitudinal velocity difference in a turbulent flow as measured on 
an intermediate length scale r. In the simplest model, we can take |m^ as a measure of 
the kinetic energy (per unit mass) of the flow associated with scale r. We, therefore, take 
X(r) = ]^v?{r) as the relevant (scale-dependent) extensive quantity (per unit mass) with its 
rate of change, X{r)^ as the (scale-dependent) current. Were intermittent fluctuations not 
present in the current, the turbulent system would relax to the equilibrium state with mean 
energy X{r) = ^{u'^{r))o, and the fully developed turbulent state would be an equilibrium 
state. The normalized variable is, thus, x(r) = u^ / {u^)q and the equilibrium p.d.f. Pq{x) = 
e~^ determines the distribution of |m| (— oo < u < oo) according to Pod'^l) = 7i?y"^ '"' ^^" ^°- 
This is the p.d.f for u irrespective of its sign, i.e. Pod'^l) = Po{u) + Pq{—u), and so the 
even part of the w-p.d.f., po{u), is only obtained; the odd part (which must satisfy the exact 



constraint imposed by the Kolmogorov's four-flfths law 2l| for the third moment) remains 
undetermined. This implies that we can obtain only the even moments of the velocity 



differences. However, the moments are sometimes defined in terms of \u\ instead of u, in 
whicfi case tlie even moments coincide. Altfiougli tliere lias been no tlieoretical justification 
for tlie spectrum of tlie odd moments tlius obtained, experimental and numerical evidence 
suggest that the spectrum does not change significantly 22|-|25l| . In particular, the third 
moment of |m| has a scaling exponent that is close to unity, which is curious since the four- 
fifths law pertains to u and not \u\. This has motivated the very useful notion of extended 
self-similarity (ESS) 22|, |2J], according to which moments are conveniently plotted against 
the third moment as a log-log plot. As shall see, using the absolute velocity differences, 
our theory drives an exponent very close to unity for the third moment, thus providing an 
explanation for the ESS. Hence (m = 1, 2, . . .) 



(InDc 



\urpo{\u\)d\u\ = r{- + l){uX^\ 
2 



(9) 



Since current fluctuations are ignored, the free parameter of the p.d.f. can be determined 



a la Kolmogorov 26|] as (u^) 



~ {^r 



l2/3 



to leading order. Then (|u|''")o ~ (er)™"/^, which is 



the expected Kolmogorov power-law spectrum. Therefore, the Kolmogorov spectrum would 
apply if the fully developed turbulent flow were an equilibrium phenomenon. 

Let us now consider adding a Gaussian white noise, which is anticipated at sufficiently 
high values of the control parameter, to the current X{r). Since the current is scale depen- 
dent, we consider scale-dependent fluctuations by taking D = D{r) so that q{r) = 1 + Z)(r)r 
now depends on the scale. In the presence of such intermittent fluctuations, the system 
relaxes in the quasi-equilibrium state of fully developed turbulence where the mean dissi- 
pation vanishes (at intermediate scales) yielding zero mean current. The quasi-equilibrium 
p.d.f. is the (scale-dependent) escort distribution ([H]), which similarly implies 



p{\u\) 



2\u\ 



1 + 



u 



q-l __ 

2-q) («2) 
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(2-g)(w2) 
It has two free parameters, namely, q and {u^)- The moments are, hence, given by 



(i«r) 



\urp{\u\)d\u\ 



1 m 1 



t) 



m/2 



/O (g — 1) 2+1 JO 

The integral exists provided m < 2{q — 1)^^, which includes sufficiently large m's when q is 
close to unity. We assume this is the case (see shortly). Then 



(|m|™) 



Vg-l 



2 ^ 



' 2-q 

r(f + 1) («2)n./2 r(-i^) V^-i, 



m/2 



(10) 



which reduces to ([9]) for g = 1, of course. This will give the velocity statistics in the fully 
developed turbulence state. Our aim is to drive the structure functions in the form 

(luD ~ (err/^iy)^- (11) 

to leading order, where L(^ r) is the integral length scale and the proportionality constants 
are non-universal, depending on the detailed geometry of the turbulent production. 
The free parameter, (u^), is to be determined by flTTl) with rra = 2, so that 

{u') ~ {erf/%^r. (12) 



From (TTOl) we have 



\un _ r(f + i).^^2^i^ ^2 - g^^ ^ r(^ - f] 



{\u\m') r(f + i)^ ' \q-ij r(^-f 

which, for m > m' both odd/even, becomes 



m — m 



Since (2 — g)^^ = Q + 0{{q — 1)^) = g to leading order, and the product is 



(13) 



exp|- E ln[l-(/ + ^)(g-l)] 
exp U^ + m) - (^ + m')]^ + 0((g - 1)^)) = g[(^+™)-(-^W)]/4 
having used Ing = (g — 1) + 0((g — 1)^), f fT3|) yields 

(|«P) = r(- + l)(M2)f g™(™-2)/8 (14) 

which checks for m = 2 and reduces to iQ for g = 1, as it should. For m = 4, we have 
(m^)/('U^)^ = 2g, so that the (scale dependent) value of the kurtosis (flatness) determines 
the other free parameter, q{r). Using (ITTl) . we find 

Qir) ~ {^)-^'^^-^^\ (15) 

Since 2'r]2—ri4 <^ 1, g(r) is a slowly varying function of r in the intermediate scales and, hence, 
so is the upper bound mmax(^) = 2(g — 1)^^ on m. Although the reliability of experimental 

8 



data rapidly decreases for moments of order higher than m = 10, a scale-dependent upper 
bound on the order is at variance with the general supposition that moments of all orders 
exist. However, as we shall see below, it will serve to remove the shortcomings of the log- 
normal model with regard to its decreasing higher order moments and the Novikov inequality 
[l9|. Substituting in (TTl]) our expressions for (m^) and q from (IT^ and (TT^ . we finally obtain 
the spectrum (ITT]) with the intermittency exponents 

Vm = --m{m - 3)r]2 + -m{m - 2){2r]2 + m)- (16) 

Z o 



The first term coincides with the prediction of the log-normal (LN) model [27|] (with 
772 = /i/9 as the free parameter), which is known to yield good agreement with experimental 
data for moments up to order m = 10. Indeed, if we substitute the LN value 774 = —2ri2 
in ( TT6|) . then the second term vanishes. Thus, the second term represents deviation from 
log-normality and yields corrections to the LN predictions. Denoting the relative deviation 
of 7/4 from its LN value by 6 {\6\ < 1), (TT6|) can be written as 

Vm = vl^ - \rnim - 2)r]25. (17) 

Writing ( 1TT1) as d'ul'") ~ r''™, the scaling exponents Cm = "^/S + rjm of our 'corrected' LN 
model are, therefore. 

Cm = C - l^m - 2)^25. (18) 

Since the LN exponents for higher moments generally fall below experimentally measured 
values, we anticipate a negative 5 (for example, 5 ~ —0.08, as given by the She-Leveque 



model 281]). For m > ttIq^ ~ I + 3^(~ 16.5) (for numerical estimates we use the LN 
value, ri2 ~ 0.2/9), the scaling exponents (^ monotonically decrease with m, implying the 
breakdown of the incompressibility approximation in the fully developed turbulent state 



29l |. From (fT8|) . Cm is a monotonically increasing function of m provided m < tuq — -^^ 



2+5 ' 
and for 5 < we have mo > mlf^ . Therefore, the exponents Cm increase with m if the slowly 

varying upper bound r?7,max('") remains sufficiently close to ttiq^. That is, if g — 1 remains 

close to 2/mQ^{^ 0.12), which is consistent with our assumption that g — 1 ^ 1. In other 

words, for \q(r) — 1 1^1 ^ 1, the corrected model (with 6 < 0) rids the LN model of 

its decreasing higher order moments, yielding structure functions only up to order [rriQ^] ([] 

stands for the integer part). Then, it will also satisfy the Novikov inequality r/m + '^ > 0, 



TABLE I: Comparison of ("„ as given by our theory, the LN model, and experiment. 



m 



Theory 



LN Model 



Ref. [24] 



Ref.'^[30] 



1 


0.36 


2 


0.69 


3 


1.00 


4 


1.29 


5 


1.56 


6 


1.81 


7 


2.04 


8 


2.24 


9 


2.43 


10 


2.59 



0.36 
0.69 

1 
1.29 
1.56 
1.8 
2.02 
2.22 
2.4 
2.56 



0.37 
0.70 
1.00 
1.28 
1.54 
1.78 
2.00 
2.23 



0.71 

1 
1.33 
1.54 

1.8 
2.06 
2.28 
2.41 
2.60 



"Moments pertain to u and so ('s = 1 by assumption. 



which is violated by the LN modeh This can be readily checked using (llTl) and noting that 

—0.08, Table [T] compares Cm up to order 



6-<5 ^ ^^ 



6m^+(5- 



10-S ^ 2 ^ 2+5 

Using the tentative values ri2 = 0.2/9 and 6 = 
m = 10, as given by our theory (equation (fTSjl ). the LN model and the references 



24 



3Q|. 



The agreement is quite good considering the uncertainties involved. In particular, we notice 
the value of Cs, which is only slightly larger than the exact value 1 (implied by the four-fifths 
law for {u^)) in conformance with {\u\^) > \{u^)\ and the data of [24|, theoretically justifying 
ESS. 
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